Skip to content

[PWGJE,Infrastructure] Create qg tree creator#15555

Open
Physicistdana wants to merge 29 commits intoAliceO2Group:masterfrom
Physicistdana:createqgTreeCreator
Open

[PWGJE,Infrastructure] Create qg tree creator#15555
Physicistdana wants to merge 29 commits intoAliceO2Group:masterfrom
Physicistdana:createqgTreeCreator

Conversation

@Physicistdana
Copy link
Copy Markdown

Jet analysis code

@github-actions github-actions bot added the pwgje label Mar 28, 2026
@github-actions github-actions bot changed the title Createqg tree creator [PWGJE] Createqg tree creator Mar 28, 2026
@Physicistdana Physicistdana marked this pull request as draft March 28, 2026 04:45
@Physicistdana Physicistdana reopened this Mar 28, 2026
@Physicistdana Physicistdana marked this pull request as ready for review March 28, 2026 21:01
@Physicistdana Physicistdana changed the title [PWGJE] Createqg tree creator [Infrastructure,PWGJE] Create qg tree creator Mar 28, 2026
@Physicistdana Physicistdana changed the title [Infrastructure,PWGJE] Create qg tree creator [PWGJE] Create qg tree creator Mar 28, 2026
@Physicistdana Physicistdana changed the title [PWGJE] Create qg tree creator [Infrastructure,PWGJE] Create qg tree creator Mar 28, 2026
@Physicistdana Physicistdana changed the title [Infrastructure,PWGJE] Create qg tree creator [PWGJE,Infrastructure] Create qg tree creator Mar 28, 2026
//------------------------------------------------
float deltaPhi(float phi1, float phi2)
{
return std::remainder(phi1 - phi2, 2.f * static_cast<float>(M_PI));
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There is RecoDecay::constrainAngle for this.

return std::remainder(phi1 - phi2, 2.f * static_cast<float>(M_PI));
}

float deltaR(float eta1, float phi1, float eta2, float phi2)
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Isn't this an existing method?

int mpdg = mom.pdgCode();

// stop at quark or gluon
if (std::abs(mpdg) == 21 || (std::abs(mpdg) >= 1 && std::abs(mpdg) <= 6)) {
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Don't use magic numbers.

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
{
return WorkflowSpec{
adaptAnalysisTask<QGTreeCreator>(cfgc, TaskName{"qg-tree-creator"})};
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why do you need TaskName?

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the review. I updated the code implementation.

This task creates a quark–gluon tagging training table from matched detector- and particle-level jets. It stores jet substructure observables, matching ΔR, and initiating-parton label. My goal is to create a compact training table for ML-based quark–gluon discrimination study

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Development

Successfully merging this pull request may close these issues.

2 participants